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Complicated cosmic string loops will fragment until they reach simple, non-intersecting ( "stable" ) 
configurations. Through extensive numerical study we characterize these attractor loop shapes 
including their length, velocity, kink, and cusp distributions. We find that an initial loop containing 
M harmonic modes will, on average, split into 3M stable loops. These stable loops are approximately 
described by the degenerate kinky loop, which is planar and rectangular, independently of the 
number of modes on the initial loop. This is confirmed by an analytic construction of a stable 
family of perturbed degenerate kinky loops. The average stable loop is also found to have a 40% 
chance of containing a cusp. We examine the properties of stable loops of different lengths and 
find only slight variation. Finally we develop a new analytic scheme to explicitly solve the string 
constraint equations. 

PACS numbers: 98.80.Cq, 11.27.+d 



I. INTRODUCTION 

Cosmic strings may be observed through their gravita- 
tional wave emission, gravitational lensing of background 
galaxies, imprints on the cosmic microwave background, 
particle emission, and various other signatures. To work 
out the signatures of cosmic strings it is essential to know 
the properties of the string network, the number density 
of loops and their length and shape distributions. Sev- 
eral observational signatures, including all burst events, 
whether gravitational or other, depend on the sudden 
whip-like motion of the string, a feature called a "cusp" . 
Yet other signatures are affected by sharp corners on a 
string, called "kinks". (See Ref. [1] for a review.) Hence, 
to make reliable observational predictions from cosmic 
strings, it is important to characterize the generic prop- 
erties of cosmic string loops and to quantify the frequency 
with which kinks and cusps occur. 

Scherrer and Press (SP) [5] undertook a study of the 
dynamics of large loops starting with randomly chosen 
loops containing 10 harmonic modes. As each loop os- 
cillates, it self-intersects and reconnects, breaking into 
smaller loops. In SP the loop is approximated by 128 dis- 
crete segments thus limiting the resolution of their study. 
The study in SP was followed up by Scherrer, Quashnock, 
Spergel and Press (SQSP) [3] with a focus on evaluating 
the gravitational power emitted from string loops. A sim- 
ilar study was undertaken by Casper and Allen (CA) [3] 
and is the most extensive study done to date. CA con- 
sidered loops with 10 harmonics and discretized the loop 
using 600 segments. 

In the present work we have built upon the SP scheme, 
improvements coming from more computational power 
and choice of algorithms. To put our current study in per- 



spective, we have considered loops with 3 to 50 harmonics 
and discretized using several times 10,000 segments. Our 
very fine resolution ensures that discretization effects are 
unimportant as will be clear from our results. A compar- 
ison of previous loop fragmentation studies with this one 
is given in Table [T] 

Our conclusions can be summarized as follows. Large 
loops with many oscillation harmonics split into smaller 
loops until the small loops have a minimal number of 
harmonics. In other words, the number of fragments is 
directly proportional to the number of harmonics on the 
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Loops 


SP 


125 


a 


10 


20 


561 


CA 


60C 
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10 


200 


5,723 


Present work 


10,000 


3 


1,000 


8,308 




10,000 


10 


3,000 


94,628 




10,000 


20 


1,000 


63,490 




10,000 


30 


1,000 


96,207 




10,000 


40 


1,000 


128,764 




10,000 


50 


1,000 


157,968 




50,000 


10 


1,000 


32,158 




50,000 


50 


1,000 


162,157 



a Loops were also run with 256 segments without producing 

significantly more small loops 
b The loops were rerun with 800 segments resulting in almost no 

new daughter loops being produced. 



TABLE I. Parameters from studies of loop fragmentation. 
Modes refers to M in Eq. Q. The numbers provided here 
are for the "Type A" loops defined by SP and used by CA; 
see the discussion after Eq. Q for more details. For SP see 



initial loop; on average 3Af loops are created from an ini- 
tial loop with M modes. The final non-self-intersecting 
( "stable" ) loops are approximately planar and their left- 
and right- moving modes tend to be orthogonal. The 
average stable loop contains four kinks and has a 40% 
chance of containing a cusp. Hence, realistic cosmic 
string loops should be visualized as oscillating rectangu- 
lar shapes, like the "degenerate kinky loops" discussed in 
Ref. |5j . This picture is strikingly clear in the animations 
available in Ref. [6]. 

The conclusion that stable loops are close to degener- 
ate kinky loops also means that gravitational and other 
radiation from stable loops can be estimated by the ra- 
diation from degenerate kinky loops, which is simple to 
calculate analytically 5J . Knowing the shape of loops can 
also help with other observational signatures e.g. lensing. 

Our detailed results can be found in Sec. lIIII At a finer 
level, the shape of a cosmic string loop may also depend 
on its length. We investigate this question in Sec. |IV| and 
find only weak dependence. It is a fair approximation to 
think of long loops as also being approximately planar 
and rectangular. In Sec. [V] we give analytical arguments 
that help us understand our numerical results. Notable is 
our demonstration of a set of perturbed degenerate kinky 
loops that are stable. We conjecture that such loops are 
attractors for the evolution of cosmic string loops. 

In the Appendices we describe technical details of the 
numerical algorithm, formulas to boost to the rest frame 
of a loop, and finally, a new scheme to explicitly solve the 
Nambu-Goto equations and the string constraint equa- 
tions. This new scheme is very suitable for constructing 
cosmic string loops but in our numerical work we use the 
SP scheme to aid comparison to previous work. 

II. PROCEDURE 

A string in flat spacetime is described by A M = 
(t,x(a,t)) where a is the parameter along the string 
and t is Minkowski time. The solution of the Nambu- 
Goto string equations yields independent left- and right- 
moving modes 

x(a,t) = -[a{a-) + b(a + )}, (1) 

where a± = a ± t. We choose equal intervals of a to 
label equal amounts of energy which gives the additional 
constraint 

\p\ = i = \q\ (2) 

where p = a' and q = b' . Closure of a loop of length L 
implies that 

pder = — I qda (3) 
Jo Jo 

while, in the center of momentum of the loop, each of 
the two integrals vanishes. A helpful geometrical picture 
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FIG. 1. When two strings intersect, they reconnect in an 
"intercommutation event with the production of four kinks 
(discontinuities in the tangent to the string), two on each 
string. 

is that the vector functions p and q trace out a path on 
the surface of a two dimensional sphere, also called the 
"Kibble-Turok (KT) sphere" . 

The velocity of a point on the string is (q — p)/2. 
Hence, at points where the curves corresponding to p 
and q intersect, the velocity of the string reaches the 
speed of light. Such momentarily light-like points on the 
string are known as "cusps". Finally note that p and q 
need not be continuous. If there is a break in either of 
the curves, it implies a discontinuity in the tangent vec- 
tor to the string, which will appear as a sharp corner or 
"kink". 

The Nambu-Goto description of a cosmic string breaks 
down when two string segments intersect. The Nambu- 
Goto dynamics has to be supplemented by the condi- 
tion that the strings intercommute, i.e. reconnect, at the 
point of intersection. This is sketched in Fig. [T] which 
also shows the four kinks that are created during inter- 
commutation, two on each string. 

To numerically model strings we follow SP and choose 
initial loops from an ensemble constructed by expanding 
the left- and right-movers in Fourier series 

M 

a i( s ) = ° m ' 1 cos ( ms + 4>a,m,i), 
M 

his) = ^2 6 m ,i cos(ms -I- (f> b , m ,i), (4) 

m— 1 

where i stands for the three components of the vectors. 
Here we have taken s £ [0, 2tt]. The mode coefficient 
vectors, a m and b m , are chosen so that each component 
of these vectors is uniformly distributed in the interval 
[0,1]. The phases, 4> a ,m.i and 4>b,m.i, are chosen to be 
uniformly distributed in the interval [0,2tt]. This corre- 
sponds to the "type A" scheme in SP and CA. Had we 
chosen the length of the mode coefficient vectors in the 
interval [0, 1/m 2 ], we would have a direct correspondence 
with the "type B" scheme of SP and CA. As already dis- 
cussed in CA, the type A loops have more power on small 
scales and leads to more fragmentations. For this reason 
we have only studied type A loops in this work. However, 
the shapes of the final stable loop population is expected 
to be the same in both schemes. We have checked this 
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by considering several different values of M since this 
controls the power on small scales. A comparison of our 
choice of parameters with those of earlier work is given 
in Table U 

The challenge now is to meet the constraints in Eq. ^ . 
This is done by recognizing that we are free to choose 
the parameter s as a function of the parameter a. By 
differentiating Eq. Q we obtain 



\P\ 



ds 
da 



M 

E 

m—1 



(5) 



This provides a differential equation for a in terms of s 
which we solve numerically with the additional require- 
ment that ds/da > 0. To insure that s(er) £ [0, 2tt] for 
a e [0, 1] we first rescale a m and b m before solving 
A numerical inversion then gives s as a function of a 
which finally allows us to obtain p(a) and q(a) at the 
initial time, t = 0. Note that the value of M corresponds 
to the "number of harmonics" on the loop in terms of 
the parameter s but not in terms of the invariant length 
a. However we will still refer to M as the number of 
harmonics. 

Once we know the function p(er) and q(a), the Nambu- 
Goto evolution is straight-forward since only the argu- 
ments of the functions change as given in Eq. (!]). The 
numerical implementation of intercommutation events is 
more involved and described in Appendix [A] 

An intercommutation causes the original loop to frag- 
ment into two smaller loops, with each of the two loops 
having one left-moving kink and one right-moving kink. 
In this way, the loops keep fragmenting with the produc- 
tion of more and more kinks. The fragmentation ceases 
when a loop reaches a non-self-intersecting configuration. 
It is these stable loops that we wish to study and char- 
acterize. 

In principle, two fragments of the initial loop can in- 
tersect and reconnect to form a bigger loop. In practice, 
since the fragments are flying apart, such collisions do not 
happen frequently. In any case, we ignore such mergers. 



III. NUMERICAL RESULTS 




FIG. 2. A sample initial loop with 10 harmonic modes. The 
darker, thicker segments of the loop are nearer the viewer. 
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FIG. 3. Length of loops versus generation. By generation 10 
the distribution is nearly identical to the stable loop length 
distribution. As fragmentation occurs the length of loops de- 
creases. The distribution of the predominantly small length 
stable loops can be seen in Fig. [4] 



Most of our data is collected for loops of 10,000 seg- 
ments and 10 harmonic modes; see Table [TJ In Fig. [2] we 
show a sample initial loop. We evolve the loop as de- 
scribed in Sec. [TT] and keep track of the fragmentations. 
We define loop "generations" in the following way. The 
initial loop is the first generation. When it splits into 
two loops, those are in the second generation, and so on 
for successive generations. Once a loop is in a stable 
configuration it is counted in each subsequent genera- 
tion. For example, a stable loop that is generated in 
the third generation will be counted in the third, fourth, 
etc. generations. At every generation, we can plot the 
length distribution of loops as shown in Fig. [3j The plot 
shows that the second generation has roughly a uniform 



distribution of lengths but that the stable loop distribu- 
tion is sharply peaked at small lengths. The stable loop 
length distributions are plotted in Fig. |3J Though very 
small loops are formed through fragmentation the peaks 
in this figure at L ~ 2/N show the resolution limit of 
the simulations. Clearly the bulk of the stable loops are 
above our resolution limit. For comparison the SP and 
CA results are shown and are in good agreement with 
our results. 

From the length distribution we calculate the average 
length per generation for several values of M. In Fig. [5] 
we show the average length scaled by M for the stable 
loops. All the plots for M = 10, . . . , 50 asymptote to the 
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FIG. 4. Length distribution of stable loops shown on a log- 
arithmic scale to better see the small length behavior. The 
resolution limit of the simulations is seen by the small peaks 
at L ~ 2/N. The length distribution is otherwise indepen- 
dent of resolution, N. Comparison to the SP (solid, gray line) 
and CA (dashed, gray line) shows good general agreement. 



FIG. 5. Average loop length multiplied by the number of har- 
monics in the parent loop versus generation. This weighted 
length asymptotes to ~ 0.33 for all initial loops. 



same value, 0.33, which clearly demonstrates that the 
average length of stable loops is proportional to 1/M. 
The plot for M = 3 asymptotes to 0.36 and suggests 
that loops with only very low harmonics could behave 
differently from loops with high harmonics. However, 
even the M = 3 value of 0.36 is very close to 0.33 and 
we will disregard this difference in the discussion below. 
Then the asymptote at 0.33 = 1/3 implies that an initial 
loop with M harmonics fragments into 3M stable loops, 
as can also be read off from Table|U Given that the initial 
loop is straight on the length scale 1/M, the stable loop 
is also composed of roughly straight segments of length ~ 
1/3M. This result solidifies and quantifies the finding in 
Ref. [7J that fragmentation does not continue indefinitely 
to smaller and smaller loops. 

In terms of the p and q curves on the KT sphere, since 
a loop with only the fundamental harmonic wraps around 
a great circle once, we expect that stable loops will wrap 
the sphere only a third of the way. Also, the discussion 
below shows that the p and q curves for the stable loops 
occur in two disconnected segments on the KT sphere. 
Hence each segment only covers a sixth of the KT sphere 
and can be thought of as an arc on the KT sphere that 
extends 60° on average. 

Our result that the stable loops are made from roughly 
straight segments can only work if the loops have kinks 
on them. It can be argued that the average number of 
kinks should be four jj |. Assume that there is an initial 
loop with Nk pre-existing kinks. (In our numerical work, 
the initial loops were smooth and so Nk = 0.) Every 
intercommutation event adds four kinks and one extra 
loop to the system. Therefore after n intercommutations 
there will be 4n + Nk kinks and n + 1 loops. So the 
average number of kinks per loop is (4n + Nk)/(n + 1) 
which goes to 4 in the large n limit. The distribution of 
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FIG. 6. The number of kinks on stable loops. Each loop must 
have at least 2 kinks. Most of the loops have between 2 and 5 
kinks independent of resolution and the number of harmonic 
modes. 



kinks, however, needs to be calculated numerically and is 
shown in Fig. [6] From this argument and the numerical 
simulations we have the picture that stable cosmic string 
loops have roughly straight sides and four corners (kinks) 
on average. 

In general, stable loops have relativistic center of mass 
velocities as shown in Fig. [7J To analyze the shape of a 
stable loop, we must first transform to its center of mass 
frame. Care must be taken to transform the coordinates 
so that the gauge conditions in Eqs. ^ are satisfied as 
described in Appendix [B] Let us denote the resulting 
curves on the KT-sphere in the rest frame by P and Q. 

To analyze the planarity of a stable loop, we calculate 



5 




0.4 0.5 0.6 
Velocity 



1.0 
0.9 
0.8 
0.7 

a) 
ot 

S 0.6 

> 

< 

| 0.5 
> 

5 0.4 
0.3 
0.2 
0.1 



,»»»«««««« 



♦ ♦ Largest Eigenvalue 

• • Middle Eigenvalue 

■ ■ Smallest Eigenvalue 



10 12 
Generation 



FIG. 7. The velocity distribution of stable loops for various 
resolutions and numbers of modes. The distribution is uni- 
versal and sharply peaked at v ~ 1. Also shown are the SP 
(solid, gray line) and CA (dashed, gray line) for comparison. 
The resolution of SP did not allow them to produce high ve- 
locity loops. Our results refine the peak at high velocity found 
by CA. 



the "moment of inertia" tensors for P and Q via 



Aij — 



B, 



J PiPjda-, 



QiQi citr 4 



(6) 



If the P (and similarly for the Q) curve is uniformly 
distributed on the KT-sphere, the three eigenvalues of 
the A tensor would be equal. If one of the eigenvalues of 
A vanishes, then P is distributed in a plane; and if two 
eigenvalues vanish, the distribution is lineal. In Fig. [8] 
we show the average eigenvalues as the loop continues to 
fragment. It is clear that one of the eigenvalues vanishes, 
and the largest eigenvalue is much larger than the middle 
eigenvalue. So the P curve is mostly in one direction 
on the KT-sphere (say around the z-axis), with a little 
spread in some other direction. Since the centroid of the 
P curve has to vanish in the center of mass frame (see 
the discussion below Eq. ([3])), the P curve must contain 
two short segments, say one near the north pole and the 
other near the south pole. In an idealized ("degenerate") 
case, the P curve would consist of just 2 points on the 
KT-sphere. Similarly the Q curve would consist of just 
2 other points. Since the loop is a sum of the left- and 
right-movers, this implies that the loop is planar and 
defined by the plane of the P and Q curves. 

As another check of this picture, we have plotted the 
kink sharpness for the p curves, defined as [8] 



-(1-P--P+), 



(7) 



where P± are the values of P on either side of the kink. 
The plot in Fig. [9] demonstrates a peak near sharpness of 



FIG. 8. The average of the eigenvalues of the P and Q "mo- 
ment of inertia" tensors Q versus loop generation. Since 
both P and Q have the same properties, both are included 
in the average to increase the statistics. We notice that the 
smaller eigenvalue approaches zero and the largest eigenvalue 
is much larger than the middle eigenvalue showing that the 
P and Q curves are localized on the KT-sphere. This plot 
is for N — 10, 000, M = 10. The asymptotic behavior is 
independent of resolution and the number of modes. 
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FIG. 9. The kink sharpness distribution (see Eq. Ff\. The 
distribution is independent of resolution and number of modes 
on the initial string. We see that the distribution is sharply 
peaked at 1, that is where the P before and after the kink are 
anti-parallel. 



1, which implies P_ sa — P + . This shows that there are 
180 degree jumps on the KT-sphere, in accordance with 
our picture that the P curve corresponds to two small 
anti-podal regions on the KT-sphere. 

We can go a bit further and ask for correlations be- 
tween the left- and right-movers. F or this we plot the 
distribution of ((P ■ Q) 2 ) in Fig. 10 (Angular brackets 
denote average over a loop.) This plot shows that P and 



Q tend to get more orthogonal with fragmentation. 
This picture of stable loops implies that the P 
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FIG. 10. The distribution of {(P ■ Q) 2 ) showing that the left- 
and right- movers become more orthogonal to each other with 
more fragmentation. 
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FIG. 11. The average number of cusps on each stable loop. 
On average we see that there is approximately a 40% chance 
of a stable loop to contain a cusp. This result is independent 
of resolution and the number of modes on the initial loop. 



— Q curves rarely cross each other and hence that cusps 
should be suppressed. For an initial loop with M modes 
we expect it to contain approximately M 2 cusps. (Each 
mode roughly corresponds to a great circle on the KT- 
sphere so the number of intersections is proportional to 
A/ 2 .) 

We have seen that on average 3M stable loops are pro- 
duced. If / is the "cusp survival fraction" — the fraction 
of cusps on the initial loop that survive on the stable 
loops — then on average we expect M 2 f/3M cusps on 
each loop. Fig. [TT] shows that, on average, each stable 
loop has a 40% of containing a cusp independent of res- 
olution and the number of modes on the initial loop so 
the cusp survival fraction is found to be / ps 6/5M. (We 
have verified this is true in our simulations but do not 
show the results here.) 



FIG. 12. The velocity of stable loops versus the length of 
the loop. The x error bars show the width of the length bins 
employed and the y error bars show the 95 percentile ranges in 
each length bin. The square symbol shows the average in the 
bin. Short loops have center of mass velocity ~ l/v2 which 
is the root-mean-squared velocity of the string in the initial 
loop. Longer loops tend to have a lower velocity consistent 
with the fact that the initial loop is at rest. 



IV. LENGTH DEPENDENT RESULTS 

The focus of this work is on the average properties of all 
stable loops. The statistics we have presented give equal 
weights to all loops. For some observational signatures, 
however, long loops may be more relevant than small 
loops. Hence, in this section, we provide some statistics 
as a function of loop length. The results shown here are 
for the N = 10, 000 and M = 10 runs. Logarithmic 
length bins are chosen based on the distribution of loop 
lengths (see Fig. [4]) . Ten bins are chosen between lengths 
of 10~ 3 and 0.1. The initial bin begins at a length of 
5 x 10 -4 to ensure that all loops are above our resolution 
limit; approximately 94% of the loops are included. The 
final bin includes all loops with lengths longer than 0.1. 

The results are shown in Figs. [l"2}fT5| These figures 
show the velocity, number of kinks, number of cusps, and 
eigenvalues, respectively, for each length bin. In the plots 
the x error bar shows the width of the bin and the y error 
bar provides the 95 percentile range for the values. The 
results are as expected. The smallest loops are predom- 
inantly high velocity, contain approximately four kinks 
per loop, have a small chance of containing a cusp, and 
are planar. The largest loops are predominantly lower ve- 
locity, contain more than four kinks, have a higher chance 
of containing a cusp, and are less planar. 



V. DISCUSSION 

Our numerical results point to the picture that sta- 
ble loops are deformations of degenerate kinky loops. In 
this section we provide further evidence for this picture 
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FIG. 13. Similar to Fig. [12] now for the number of kinks 
on stable loops versus the length of the loop. Longer loops 
tend to have a larger number of kinks though the variation 
is slight. The loop with the maximum number of kinks has 8 
kinks instead of the canonical 4. 




FIG. 14. Similar to Fig. [12] now for the average number of 
cusps on stable loops versus length of the loop. Longer loops 
typically contain more cusps on average. The loop with the 
maximum number of cusps has only 3 cusps. 



by demonstrating that there exists a class of perturbed 
degenerate kinky loops that are stable. 

We begin by describing the degenerate kinky loop [5], 



Pdk = cos(a )z, 
<7dk = cos(/3 )x, 



where 



a = 7rL2er_J. 



(8) 



(9) 



[x\ denotes the greatest integer less than or equal to x, 
and a± <E [0, 1]. In Fig. 16 we show some snapshots of a 
degenerate kinky loop. 

Note that the degenerate kinky loop collapses to a dou- 
ble line twice during an oscillation period, yet the stable 



FIG. 15. Similar to Fig. [12] now for the eigenvalues of the 
moment of inertia tensors of the stable loops versus length. 
Longer loops are also planar and only slightly less so than 
smaller loops. 




FIG. 16. Snapshots during the evolution of the degenerate 
kinky loop. Note that the loop periodically collapses to a 
double line. 



loops do not self- intersect. This must be due to the fact 
that stable loops are not exactly degenerate. To show 
that this is the case, consider 

p = cos(a)z, 

q = cos(27T6er + + p)x + sin(27reer + + /3)y (10) 
where < e < 1 and now 

a = ir[2a-\ , 

/3 = (l-eM2-7 + j. (11) 

As illustrated in Fig. [TTJ these perturbations stretch out 
the point corresponding to q on the KT-sphere for the 
degenerate kinky loop to an arc of length en [5] . 

We can check explicitly that the self-intersection con- 
dition 



x(a, t) = x(a' , t) 



(12) 



is satisfied for such a loop if and only if a = a' (mod 
1), for any t, that is, the loop does not intersect at any 
time. We first write the intersection condition in terms 
of a and b, 



i(cr -t) + b{a + t) = a(a -t) + b(a' + t). 



(13) 



Since for our loop \10\ we have p ■ q = for all a and t, 
we also have a ■ b = 0, and self-intersection requires 



a(a — t) = a((i' — t), 
b(a + t) = b(a' + t). 



(14) 
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FIG. 17. The p and q curves for a perturbed degenerate 
kinky loop ( |10[ ) plotted on the KT-sphere. The p curve is 
represented by the two points at the poles. The q curve is 
represented by the arcs along the equator. 



The first condition is satisfied for er_ = a'_ or for a'_ = 
1 — cr_ . Only the latter possibility is of interest since we 
are looking for an intersection with a =/= a' . The second 
condition for the x and y components leads to 

sin(2e7rcr+) = sin(e7r) — sin(2e7r<7^_ — en) 
cos(2e7T(T + ) = 1 + cos(e7r) — cos(2e7ro/|_ — en) (15) 

Note that both equations must be satisfied for a self- 
intersection to occur. 

With the help of some trigonometric identities, we can 
rewrite (fl5l) as 



cos 



C7T I (7 _i_ — (7 , 



(?)«»(?)• 



en I rr 

,2 



(?) 



(16) 



Taking the ratio of these equations leads to 



tan 



tan 



so that o' + + <t_|_ = 1. By squaring and adding the equa- 
tions, we obtain 



cos 



± cos ( — 



<(?) 



(18) 



so that the second relation in Eq. ( 16 1 leads to 



cos 



cos 



£7T (7+ 



or 



sin (e7r(T+) sin 



C7T 0" , - 



= 



(19) 



(20) 



which implies cr+ = or else a' + = 1/2. (Solutions such 
as er + = 1 je can be ignored because we restrict <r + to lie 
in the interval [0, 2tt] and |e| < 1.) Therefore the only 



solution is of the type a + = 0, <r' + = 1 or else a' + = 1/2, 
<t+ = 1/2, which are both trivial. This shows that there 
are no self-intersections. The conclusion is intuitive since 
b is obtained by integrating q' and hence is a vector that 
lies in the xy— plane (see Fig. 17) and rotates in this 
plane. Since < e < 1, b never attains the same vector 
value twice for a + e [0,1]. 

The perturbed degenerate kinky loop discussed above 
does not contain any cusps while our numerical results 
show that roughly half of the stable loops have a cusp. 
To understand the P and Q curves of a stable loop with 
cusps we show an example of these curves for a sample 



loop from our numerical runs in Fig. 18 which contains 
three kinks and one cusp. In this example the P curve 
is like our perturbed degenerate kinky loop but the Q 
curve is much more elongated. The two curves are in 
planes nearly perpendicular to each other. 

To construct a stable loop with cusps consider 

p = sin(/3_)x + cos(/3_)z, (21) 
q = sin(/3 + ) [sin ( <f> + ) x + cos(<fi+)y] — cos(f3 + )z, 

where 



0± =2ire ± <j ± + (l-e ± )irl2a ± \, 



(22) 



and <p + is an arbitrary phase. There are two cusps due 
to the intersections at ±z. These occur at a = t = and 
a = t = 1/2. We have numerically confirmed that this 
loop is stable for the parameter choice 



1 



1 



58' 



95' 



67 



(23) 



In Fig. [18] we show the p and q curves for a stable loop 
with cusps similar to (21 ) extracted from our simulations. 



VI. CONCLUSIONS 

We have shown that the fragmentation of a large 
complicated loop yields stable quasi-planar, quasi- 
rectangular loops, that are similar to degenerate kinky 
loops. If cosmic string loops are the result of a large 
number of intercommutings - even as few as five or six 
intercommutings may be sufficient (e.g. Fig. [8]) - we ex- 
pect that they too will be quasi-rectangular. Earlier work 
has examined the stability of relatively simple loops (with 
a small number of harmonics) and found stability over a 
region of parameter space [9l412j . It is unclear if cosmo- 
logical evolution can directly produce these stable loops 
with few harmonics. Cosmological loops will certainly 
contain kinks and this feature is missing in the loops in- 
vestigated in these earlier studies. 

In a cosmological setting, when a loop is produced 
from the infinite string network, it will fragment down 
to stable loops within two oscillation periods, and there- 
fore only stable loops are relevant for cosmological signa- 
tures. Indeed the aim of the CA and SQSP studies was 
to obtain the gravitational power emitted from realistic 
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f. • • • • 



FIG. 18. The P and -Q curves plotted on the KT-sphere for 
a stable loop from our simulations. This loop contains three 
kinks and one cusp. The two curves are represented by the 
dots. The Q curve is the dots along the equator. The P curve 
is the dots perpendicular to the equator (running vertically 
on the left of the sphere) and includes the lone anti-podal 
point (on the right of the sphere). Note that the points on 
the sphere are not weighted equally and the centroids of both 
P and Q curves are at the center of the sphere. 



cosmic string loops. We can also estimate the gravita- 
tional power based on the analytical results of [5] for 
emission from degenerate, kinky loops, and obtain an es- 
timate that is within 13% of the CA result. The gravita- 
tional power radiated from degenerate kinky loops is an- 
alytically calculated to be 64 ln(2) G/j 2 s» 44G/i 2 whereas 
CA numerically estimate 39G/1 2 from their simulations, 
where fi is the string tension. 

For cosmic strings, Hubble expansion, frictional damp- 
ing and radiation damping also come into play, though 
the effects are generally on very long time scales and de- 
pend on the environment in which the strings are placed. 
Hence we expect that when a large loop is produced from 
the network, it will fragment into stable loops within two 
oscillation periods and then, as these stable loops oscil- 
late for many oscillation periods, damping effects start to 
play a role and to change the Nambu-Goto dynamics. So 
to characterize the effect of damping on the cosmic string 
network loops, we can limit our study to the damping ef- 
fects on the perturbed degenerate kinky loops and the 
stable loops with cusps described above. 
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FIG. 19. Geometric splitting of a loop. Shown in this figure 
are the a± grids for various loops. The diagonal lines are 
lines of constant cr_ (up and to the right) and constant a+ 
(up and to the left). Time runs upward and a rightward. 
The left and right edges of the grids are identified to describe 
a loop of string. The initial loop is shown on the top with 
the intersection occurring at points a and b. The dashed 
lines show the constant values of a± at which the intersection 
points occur and so at which kinks are formed. The daughter 
loops are formed by cutting the grid and identifying the new 
edges. In the figure the gray shaded region is cut from the grid 
and becomes one loop whereas the remaining unshaded pieces 
are connected. The new loops are shown below the initial 
loop. The past light cone of the intersection point has been 
removed since it no longer contains useful information. The 
two intersecting diamonds are replaced by the three shaded 
diamonds in the figure. 



stitute for Advanced Study. The numerical simulations 
were performed on the facilities provided by the Case ITS 
High Performance Computing Cluster. 



Appendix A: Loop Intercommutation 

The numerical evolution in this work was performed 
using the "diamond method" for string evolution which 
is discussed in |13j . Full details of our implementation 
along with source code and documentation is available 
online [B]. Here we will discuss one aspect of the evolu- 
tion: loop intercommutation. 

The diamond method gets its name from the a± grid 
shown in Fig. [19] The top of the figure shows the grid 
structure for an initial loop. Our initial loop segments 
are of uniform length, Act, producing a regular grid. 
This grid provides a simple geometric picture of the loop. 
With this grid structure, a starting point, x(0,Q), and 
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values for p(a_) and g(cr+), taken to be piecewise con- 
stant functions, we can calculate the coordinate of the 
string x(<7, t) during its evolution if intersections do not 
occur. 



In Fig. 19 we consider the case of an intersection at 
points a and b, that is, x(a a ,t- m t) — x(ab,t- m t)- In 
this picture the intercommutation and production of two 
loops is geometric. The shaded area is excised from the 
parent loop to become one daughter loop (bottom right 
grid) and the remaining portions of the parent loop are 
reconnected (bottom left grid) . The past light cone of the 
intersection point no longer contains meaningful informa- 
tion and has been removed from the new grids. When 
an intersection occurs the grid structure is modified; an 
originally uniform grid becomes non-uniform, but other- 
wise contains the same information. In the bottom left 
grid the light gray shaded diamonds are the remnants 
of the original diamonds that intersected. The dark gray 
shaded diamond is a new one created by the reconnection 
of the grid after the middle portion has been removed. 
These two grids now describe independent loops which 
may be evolved in the usual way for all times t > t int . 



Appendix B: Loop Rest Frame Transformation 

The time coordinate of a loop, t, was chosen to be the 
proper time and the same as the background spacetime, 
a; = t; see Sec. [n] To make p and q easier to boost we 
construct four-vectors from them with the time compo- 
nents give by p° — — 1 and q° = 1. At a fixed time, t, we 
may then write in the initial loop's rest frame 



P' 



da^ 
daZ 



P 



(Bl) 



where \p\ = 1. 

The center of mass velocity of a loop is given by 



q da + — I pda 



1 



= / qd<7 + = — \ pda 



(B2) 



With 7 = 1/y 1 — boosting to the rest frame of this 
loop we find 



P>_) = 



-7(1 + v cm ■ p) 
p + (7 - l)(v cm ■ p)v cm /v* m + <yv c 



Notice that the time component of this four-vector is 
no longer —1. To correct this we apply a gauge transfor- 
mation to the coordinate er_. Let 



so that 



CT_ = 7(er_ + v cm ■ a), 
da. 



da. 



j(l + v cm -p). 



(B4) 
(B5) 



An integration gives <r_ € [0,1/7] and a coordinate 
rescaling can be used to bring the interval back to [0, 1]. 
With this we now have 



pV_) 



dd° 
da. 



d~a° 
daZ 



pu 



7(1 



■P) 



da. 
da~. 

= -1, 



(B6) 



where we have used (B3) for P. Applying the same 



transformation to the spatial piece of (|B3|) we find 

'71>cm 



P(<7_) 



p+ (7- l)(« cm 'P)Vcnx/v, 



2 

cm 



7(1 



■P) 



(B7) 



Again this looks messy but we can verify that \P\ 
as required. We can further verify that 



P{a^)da_ = I P(a_ 



da- 
daZ 



da. = 0. 



(B8) 



For q we proceed in the same way. In this case the 
gauge transformation is <r + = 7(0 + — v cm ■ b) and we find 



Q(^ 



9 + (7 - l)K:m • q)Vcm/Vcm ~ 1 V - 

tC 1 - '"cm • q) 



(B9) 



Once again we can show that \Q\ — 1. 

The above equations are not in the best form for nu- 
merical evaluation. If |t> cm | s» 1 then 7 ^S> 1 so large 
number will be subtracted from each other in the numer- 
ator. To correct this it is better to write the equations 
in terms of 7 _1 = y/ 1 — v^ m 
provides the alternative forms 

„ 7" 1 P+(1-7" 1 )K: 



since < 7^ < 1. This 



' pHm/^m + V, 



1 



P 



Q = 



7 1 Q + (l-7 1 )(w C m • q)v c m/v, 



2 

cm 



-(BIO) 



Appendix C: Explicit Solutions of the String 
Equations and Constraints 

Here we describe an explicit analytical solution of the 
Nambu-Goto equations of motion, Eq. ([!]) , and the string 
constraints, Eqs. (|3). By choosing the decomposi- 
tion in terms of left- and right- movers we can solve the 
Nambu-Goto equations as in Eq. 0. To solve the con- 
straint in Eq. ^ we start with 

p(a-) = (sin 6 cos 4>,sm 9 sine/), cos 9) (CI) 

where 9 and <j> are functions of c_. We similarly choose 
q(a+) with independent angular functions. 

We require that p be periodic under <t_ — > cr_ + 1. 
(Throughout this section we consider a loop of length 1.) 
Therefore we write 

00 

0(cr_) = 27r7_(7_ + V" 9 m cos(27rTO<j_ + a m ) (C2) 
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where j- is an integer, and 9 m and a m are arbitrary con- 
stants that can be chosen randomly to generate random 
loops. Similarly, 

oo 

</>(ct_) = 27rfc_(T_ + \_. $m cos(27rmer_ + f3 m ) (C3) 



where fc_ is an integer, and <f> m and /? m are arbitrary 
constants. 

Next we come to the closure condition in Eq. ^ . The 
integral of p will, in general, not vanish. To correct this, 
define 



J p((r_)der_ 



(C4) 



and boost p to velocity t>_ followed by a gauge trans- 
formation as described in Appendix [B] This gives us the 
final solution, 

P(^) = P+(7 -' 1)( "-' P) "- / "- +7 -"- (C5) 
V ; 7_(l + v_-p) V ; 



where 7_ = (1 — v^l) and 



d<7_ = 7_ (1 + t>_ • p)d<7- 



(C6) 



As in Appendix [5] an integration shows that er_ G 
[0, 1/7] and a coordinate rescaling can be used to bring 
the interval back to [0, 1]. 

It is straightforward to check directly that \P\ = 1 and 
also 



J P(a_)da_ = 



(C7) 



In a similar way we construct 

This scheme has the advantage that it solves the con- 
straints exactly and explicitly though it still requires inte- 
grating p(<J-) to find V- and a^. We have not employed 
this scheme in the work reported here but have remained 
with the SP algorithm to facilitate comparison. 
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